Two-site entropy and quantum phase transitions in low-dimensional models 
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We propose a new approach to study quantum phase transitions in low-dimensional lattice models. 
It is based on studying the von Neumann entropy of two neighboring central sites in a long chain. 
It is demonstrated that the procedure works equally well for fermionic and spin models, and the 
two-site entropy is a better indicator of quantum phase transition than calculating gaps, order 
parameters or the single-site entropy. The method is especially convenient when the density-matrix 
renormalization-group (DMRG) algorithm is used. 
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The search for ground state and the study of quan- 
tum phase transitions (QPTs) is a challenging problem 
when strongly correlated fermionic or spin systems are 
considered. Since exactly solvable models are rare, in 
most cases the relevant part of the excitation spectrum, 
the order parameters characterizing the various phases, 
or eventually susceptibilities are determined numerically 
on finite chains and their thermodynamic limit is deter- 
mined using the standard finite-size scaling method. Un- 
fortunately in several cases no definite conclusions can be 
drawn even if the calculations are performed on rather 
long chains. 

In this letter, we propose a new approach to detect 
QPTs and to locate the quantum critical point in low- 
dimensional spin or fermionic models. It is based on 
studying the behavior of the von Neumann entropy of 
two neighboring sites in the middle of a long chain, which 
can be defined both for fermionic and spin models, and 
can be especially easily implemented when the density- 
matrix renormalization-group (DMRG) algorithm Q] is 
used. 

The method is closely related to concepts in quantum 
information theory, which recently have attracted great 
attention in relation to QPTs. Wu et al. 3 have shown 
that quite generally QPTs are signalled by a discontinuity 
in some measure of entanglement in the quantum system. 
One such measure is the concurrence Jiwhich has been 
used by a number of authors 0, 0, 00 E3 EtJ in their 
study of spin models. Since the concurrence is defined 
for spin-1/2 systems only, for higher spins or fermionic 
models another measure of entanglement is needed. 

The local measure of entanglement, the one-site en- 
tropy, which is obtained from the reduced density matrix 
Pi at site i, has been proposed by Zanardi 11] and Gu et 
al. 01 to identify QPTs. Contrary to their expectation, 
in many cases, this quantity turns out to be insensitive 
to QPT. As an example let us consider the most general 
isotropic spin-1 chain model described by the Hamilto- 
nian 
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(1) 



In ID, the model can be solved exactly at 9 = ±7r/4 and 
9 = ±37r/4, and is known to have at least four differ- 
ent phases 13] . The ground state is ferromagnetic for 
9 < — 37r/4 and 9 > tt/2, while in between the integrable 
points separate the Haldane phase, which exists in the 
range — 7r/4 < 9 < 7t/4, from the dimerized and the quan- 
tum spin nematic phases, respectively. The existence of 
another phase, the quantum quadrupolar phase [14j near 
9 = — "in /A is not settled yet [l5(. These phases and the 
corresponding QPTs do not show up at all in the local 
entropy. As can be seen in Fig. ^ the site entropy has 
the value s, = In 3 in the whole range — 37r/4 < 9 < tt/2. 
This is due to the unbroken SU(2) symmetry in the non- 
polarized states. In the totally aligned ferromagnetic 
state, finite value is obtained for the single-site entropy, 
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1: 9 dependence of the single-site entropy of the cen- 
tral site of the isotropic spin-1 chain, determined for different 
chain lengths. 

The aim of this letter is to point out that a better 
indicator is obtained if instead of the entropy of a single 
site, the entropy of two neighboring central sites, 



Si,i+i = -Trpi,i+i In Pj. 



i+i 



(2) 



is studied, where Pi^+\ is the reduced density matrix of 
the two sites. A related quantity has been used in 0] 
to characterize the interaction between states in quan- 
tum chemistry. We will demonstrate that if — in order to 
avoid boundary effects — sites in the middle of long chains 
are considered, i.e., for i = N/2 and i = N/2 + 1, where 
N is the number of sites in the chain, anomalies are de- 
veloping in this quantity, and their analysis can be used 



2 



to detect quantum phase transitions in low-dimensional 
lattice models. This procedure is especially convenient 
when a dimerization transition is studied. In this case, 
namely, the difference of two-site entropies on neighbor- 
ing sites in the center of the chain, 



D s (i) — Sj+l,i+2 — SiA+l 



N/2 



(3) 



is finite, but the method works for other types of transi- 
tions, as well. 

As usual in the DMRG approach, we considered open 
chains. The numerical calculations were performed on 
finite chains up to 400 lattice sites using the dynamic 
block-state selection (DBSS) approach [17| . The thresh- 
old value of the quantum information loss \ has been set 
to 10~ 8 . Similar anomalies appear in the two-site en- 
tropy when periodic boundary condition is used, except 
that D s is meaningless in this case. 

As a first example the results obtained for the biliniear- 
biquadratic model are shown in Fig. [21 The upper panel 
shows SiA+i for i = N/2 and i = N/2 + 1, while D s (N/2) 
is displayed in the lower panel. 





FIG. 2: (a) 9 dependence of the two-site entropy s^i+i for 
i — N/2 and i — N/2 + 1 of the isotropic spin-1 chain for 
different chain lengths. The results for i — N/2 are displayed 
for N = 400 only, (b) Dimerization of the two-site entropy. 

In the range — 37r/4 < 8 < tt/2, one can clearly dis- 
tinguish three regimes. The most pronounced feature 
is that the two-site entropy is strongly dimerized for 
— 37r/4 < 8 < — t/4, which reflects the dimerization of 
the valence-bond structure in this regime. The quantity 
8n+i for i = N/2 and i — N/2 + 1 measures the en- 
tropy of a pair of sites coupled by strong or weak bonds, 



respectively. The lower panel shows that as longer and 
longer chains are considered, the boundaries of the dimer- 
ized region scale to the integrable points, 8 = — 7r/4 and 
-3tt/4. 

In the Haldane phase, where all neighbors are coupled 
predominantly by one valence bond, as in the AKLT state 
jl8j, the two-site entropy is not alternating, its minimum 
is at the AKLT point. Right at 8 — 7r/4, where the model 
has SU(3) symmetry, the two-site entropy is discontin- 
uous, moreover in the thermodynamic limit, a kink ap- 
pears between the left and right sides of this point, which 
is a clear indication of the transition from the Haldane 
phase into the quantum spin-nematic phase. 

Having demonstrated the usefulness of the study of the 
two-site entropy on a model where the quantum critical 
points were known, we now use this procedure to study 
the phase diagram of two ID fermionic models proposed 
to explain the neutral-ionic (N-I) transition first observed 
in organic mixed-stack char ge-t ransfer salts [lflj . 

The ionic Hubbard model |2j| is defined by the Hamil- 
tonian 



A 



(4) 



When the number of electrons is exactly equal to the 
number of sites, the competition between the on-site en- 
ergy difference A and the Coulomb energy U will deter- 
mine whether the system is a band insulator (U < A) 
where both the charge and spin gaps are finite or a cor- 
related Mott insulator (U > A) where only the charge 
gap is finite. When hopping is negligibly weak, the tran- 
sition takes place at U — A. As was first pointed out by 
Fabrizio et al. 21], for finite t values the transition be- 
tween these two states occurs in two steps. The charge 
gap closes and reopens at a critical value U c , while the 
spin gap vanishes at a different value U s . The transition 
at U = U c is Ising-like while the one at U s is a Kosterlitz- 
Thoulcss transition. A dimerized bond-order phase exists 
between the two critical values U c < U < U s . 

Since then the model has been studied in detail by 
several groups [H [H HHH H3 using both ana- 
lytic and numerical methods. Despite the many efforts 
the controversy concerning the number of transitions was 
not resolved until recently [H |H Ull The difficulty 
can be seen from Fig. |31 which shows the dimer order 
parameter (the energy difference of neighboring bonds) 
for finite chains with open boundary condition. Using 
the standard finite-size scaling procedure the vanishing 
of the dimer order cannot be established with certainty 
from these data. 

Moreover, in contrast to what has been found in the 
extended Hubbard model Q, the single-site entropy is a 
rather smooth curve without sharp extremum or discon- 
tinuity in the derivative. If, however, the two-site entropy 
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FIG. 3: (a) Dimer order parameter of the ionic Hubbard 
model for t/A — 1 for different chain lengths. (6) Entropy 
of the central site. 



is analysed, two transitions can be identified. The results 
for t/A — 1 are shown in Fig. 0] The two curves for fixed 
N correspond to s^j+i with i = N/2 and i = N/2 + 1. 
Similar behavior has been found for other values of t/A, 
too. 




FIG. 4: Two-site entropy s M+i for i = N/2 and i = N/2 + 1 
of the ionic Hubbard model for t/A — 1 for different chain 
lengths. 

The Ising-like transition, the point where the two-site 
entropy becomes dimerized, can be rather well located. 
On the other hand the Kosterlitz-Thouless transition at 
larger values of the Coulomb coupling, where the dimer- 
ization disappears again, cannot be established with cer- 
tainty by looking at D s alone. Standard finite-size scaling 
does not give a vanishing D s . However, the maxima in 
s^i+i for i = N/2 + 1 develop into a cusp in the ther- 
modynamic limit, and this should be attributed to the 
second transition. 

As a next example we have considered another model 



proposed to describe the neutral- ionic transition |2jj . In 
this model donor molecules with ionization energy / and 
acceptor molecules with electron affinity A, respectively 
alternate. If the energies are measured with respect to 
the neutral state, in which the donors are doubly occu- 
pied and the acceptors are empty, the molecular levels 
are described by the Hamiltonian 



Hq 



I V(2 - mw) + [/ D V(1 - na,yr)(l - 



A y^n 2 j-i + Ua / J"2i-i,T n 2*-i,J, ) ( 5 ) 



and the hopping, the charge transfer between neighboring 
donors and acceptors is given by 



(6) 



In realistic charge-transfer salts the on-site Coulomb 
energy is presumably the largest energy, so it is reason- 
able to assume that its unique role is to forbid doubly 
ionized (empty) donors and doubly ionized (doubly oc- 
cupied) acceptors. The transition between the neutral 
and ionic phases is driven by the intersite Coulomb cou- 
pling, 



H c = -V^(2 - n 2 i){n 2 i-i + n 2 ,+i) . 



(7) 



The relevant parameter is V/A, where A = h(I — A). 

This problem has been studied by exact diagonaliza- 
tion on relatively short chains |2^| as well as by valence- 
bond techniques (2^|. For small transfer integral the 
model was found to exhibit a first-order phase transition 
form a neutral to an ionic phase as a function of V/A. 
For t/A larger than 0.3 the transition changes character 
and a second-order phase transition has been reported. 
Our numerical result for the site entropy for t/A = 1 is 
shown in Fig. |3J A cusp is clearly seen in the one-site 
entropy at V/A ~ 1.42, indicating a single second-order 
transition, in agreement with earlier expectations. 




1.4 1.6 

V/A 

FIG. 5: Single-site entropy of the donor-acceptor model for 
t/A = 1 for different chain lengths. 

Contrary to this the two-site entropy shown in Fig. 
exhibits the same behavior as the ionic Hubbard model. 
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In fact, as will be shown in a separate paper, the two 
models are limiting cases of a unified model of neutral- 
ionic transition. Although, as for the ionic Hubbard 
model, it is not easy to locate the vanishing of dimcr- 
ization when D s alone is considered, the peak of s^j+i 
for i = N/2 + 1 produces a cusp in the thermodynamic 
limit, indicating a second transition. For small transfer 
integral, t/A < 0.3, the two transition points coalesce, 
both the one and two-site entropies exhibit a sharp jump, 
the transition is of first order. 



al. [3l|. Our experience shows, however, that calculating 
the block entropy with the same accuracy as that of two 
sites is a numerically much more demanding task. 




FIG. 6: Two-site entropy of the donor-acceptor model for 
t/A = 1 for different chain lengths. 

In conclusion, we have shown that the entropy of two 
central sites in a long open spin or fermionic chain, 
which can be easily calculated in the density-matrix 
renormalization-group procedure provides us with extra 
information not contained in the single-site entropy. The 
two-site entropy displays maxima or break points at QPT 
even in such cases when it is difficult to establish the 
opening or closing of a gap, or the vanishing of an order 
parameter, or when the single-site entropy is featureless 
since it is insensitive to the breaking of symmetry that 
distinguishes the two phases. Using this procedure we 
have demonstrated that the ionic Hubbard model and 
the donor-acceptor model do in fact have two phase tran- 
sitions. 

The procedure can be easily extended to considering 
the entropy of not just two-sites, but that of a longer seg- 
ment, whose length dependence is different whether the 
model is critical or not 30] . Fig. [7| shows the entropy of 
the block of N/2 and N/2+l sites of an N = 200 chain for 
the bilinear-biquadratic S = 1 model. The dimerization 
appears in the block entropy, as well, and the anomalies 
whose location scales to the transition points 6 — ±7r/4 
are more pronounced than in the two-site entropy. We 
have learned, after the submission of our paper, that the 
study of the block-block entanglement in connection with 
quantum phase transition has been proposed by Deng et 




FIG. 7: Entropy of blocks of N/2 and N/2 + 1 sites of the 
bilinear-biquadratic S = 1 model for a chain with N = 200 
sites. 
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